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We review the application of random matrix theory (RMT) to chiral symmetry in QCD. Starting from the 
general philosophy of RMT we introduce a chiral random matrix model with the global symmetries of QCD. 
Exact results are obtained for universal properties of the Dirac spectrum: i) finite volume corrections to valence 
quark mass dependence of the chiral condensate, and ii) microscopic fluctuations of Dirac spectra. Comparisons 
with lattice QCD simulations are made. Most notably, the variance of the number of levels in an interval containing 
n levels on average is suppressed by a factor (log n)/-K 2 n. An extension of the random matrix model model to 
nonzero temperatures and chemical potential provides us with a schematic model of the chiral phase transition. 
In particular, this elucidates the nature of the quenched approximation at nonzero chemical potential. 



1. INTRODUCTION 

Random matrix theories have been applied 
to many areas of physics ranging from nuclear 
physics |lj to quantum gravity ^ and neural 
networks ||. In general, one can divide RMT 
applications into two different groups. First, as 
a description of universal fluctuations of an ob- 
servable expressed in terms of its average value. 
For example, the Hauser-Feshbach formula ||,^| 
and universal conductance fluctuations Sec- 
ond, RMT can be used as a schematic model for 
problems involving disorder. Well-known exam- 
ples in this category are the Anderson Model , 
the Gross- Witten model || and models for neural 
networks ||]. 

In this lecture we will review the application 
of RMT to chiral symmetry in QCD. Before dis- 
cussing both types of applications, we introduce 
the concept of universality within this context. 

As has been argued, in particular by Leutwyler 
and Smilga || , the mass dependence of the QCD 
partition function in the 'mesoscopic' range 



i « V 1 / 4 « -^L, 



(1) 



where m is the quark mass, and A is a typical 
hadronic scale, is given by 

„Rc VS Tr MUe i6 ^ N f /g) 

Here, M is the mass matrix, 9 is the vacuum an- 



Z eS (M,9) = I dUe* 

'UGSU(Nf) 



gle, and S = \(qq)\- However, QCD is not the 
only theory that can be mapped onto this effec- 
tive partition function. In section 3, we will in- 
troduce a random matrix model that can be re- 
duced to (||) as well. This allows us to formulate 
universality: different theories with the same low- 
energy effective partition function have common 
properties. 

2. DIRAC SPECTRUM 

Our main focus will be on the spectrum of the 
Dirac operator, 



ijD<fr k = \ k <t>k, 



(3) 



with spectral density p(X) = J2S(X — Afe). The 
density of small eigenvalues is related to the chi- 
ral condensate by means of the Banks-Casher for- 
mula M 



1(^)1 = 



7TP(0) 
V 



(4) 



where it is understood that the thermodynamic 
limit is taken before the chiral limit. Since 
{75,^7-D} = 0, the nonzero eigenvalues occur in 
pairs ±Afc. The smallest eigenvalue is of order 
Amin ~ n/TiV. In the QCD partition function, 
the Dirac eigenvalues and the quark mass occur 
only in the combination m 2 + X 2 ,. Therefore, it 
is natural to expect that universal features of the 
Dirac spectrum can only be found in a region con- 
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sistent with (|l]) which, in terms of the eigenvalues, 
can be expressed as |A&| <C \AmmA <C A. 

In order to identify a universal quantity, let us 
consider the mass dependence of the QCD parti- 
tion function in the range (Q), 



n*./(A2 



(n fe ,/A2) 



Z eff (M, 9). 



(5) 



Here, the average (• • ■) is over gauge field configu- 
rations weighted by the QCD action. By expand- 
ing both sides in powers of we find an infinite 
family of sum- rules H . The simplest sum rule, 
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has been verified for the instanton liquid model 
of the QCD vacuum jl]J . If we write the sum in 
(g) as an integral over the spectral density and 
introduce the microscopic variable u = AUE, we 
find 



1 I u ^ du 



1 



(7) 

U£ r WE J u 2 ANf ' K ' 

What enters in (0) is the microscopic spectral 
density, which in the range A <C %/ A m i n A can be 
approximated by its thermodynamic limit 



Psiu) = lim p( — — ^ 



(8) 



The existence of this limit follows form the spac- 
ing of the eigenvalues near zero virtuality as AA ~ 
1/V in the broken phase. 

Our conjecture is that ps(u) is a universal func- 
tion that can be obtained from a random ma- 
trix model with the global symmetries of QCD. 
Note that ps(u) is not fixed by Z e g(M, 9). Before 
defining this model we discuss the well-known ap- 
plication of RMT to spectra of complex systems. 



3. SPECTRAL CORRELATIONS 
COMPLEX SYSTEMS 



OF 



RMT has a long history of successes in explain- 
ing the statistical properties of nuclear spectra 0] . 
More recently, spectral correlations have been in- 
vestigated in the context of quantum chaos |l2) . 
The starting point is the observation that the 



scale of variations of the average spectral den- 
sity and the fluctuations of the spectral density 
separate. This allows us to unfold the spec- 
trum. From the original spectrum {A^} we con- 
struct an unfolded spectrum {A' fc } with average 
spectral density equal to 1. This is achieved by 
J Xk p(X)d\ — Afe, where p(X) is the average spec- 
tral density. 

A variety of statistics has been introduced to 
analyze the spectral correlations of the unfolded 
spectrum. The best known statistic is the nearest 
neighbor spacing distribution P{S). A theoreti- 
cally simpler statistic is the distribution of the 
number of levels, nfe(n), in an interval of length 
n. In this lecture, we will consider its variance 
denoted by £2(^)5 and the As(n) statistic |l3| 
obtained by a smoothening of £2(71). 

The above statistics can be obtained analyt- 
ically for the invariant random matrix ensem- 
bles. They are ensembles of Hermitean matrices 
with probability distribution given by P{H) ~ 
exp -TiWH. Three different invariant random 
matrix ensembles can be constructed: the Gaus- 
sian Orthogonal Ensemble (GOE) when H is real, 
the Gaussian Unitary Ensemble (GUE) when H is 
complex and the Gaussian Symplectic Ensemble 
(GSE) when H is quaternion real. They are char- 
acterized by the Dyson index /3, which is equal to 
1, 2 and 4, respectively. The relevant random ma- 
trix ensemble is determined by the anti-unitary 
symmetry of the Hamiltonian H . If there is no 
anti-unitary symmetry the spectral correlations 
are given by the GUE. If the anti-unitary sym- 
metry operator A satisfies A 2 — 1, it is possible 
to find a basis in which H is real, if A 2 = —1, one 
can construct a basis in which H becomes quater- 
nion real Jl|. Spectral correlations are given by 
the GOE and the GSE, respectively. 

For simplicity we only give approximate results 
for the above statistics. The nearest neighbor 
spacing distribution is well approximated by the 
Wigncr surmise given by P(S) ~ S" ex-p(—apS 2 ), 
where ap is a constant. For n > 1 the number 
variance is given by £2(71) ~ (2/ir 2 f3) logn and 
the A 3 -statistic A 3 (n) ~ S 2 (n)/2. The RMT re- 
sults for these statistics should be contrasted with 
the results for uncorrelated eigenvalues. Then, 
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P(S) = S, £ 2 (n) = n and A 3 (n) = n/15. In 
particular, we wish to emphasize that long range 
spectral correlations are strongly suppressed with 
respect to uncorrelated eigenvalues. 

Spectral statistics of a large variety of systems 
have been compared to RMT. We wish to men- 
tion, the nuclear data ensemble [0, the zeros 
of Riemann's ^-function |Q , the correlations of 
resonances in electromagnetic resonance cavities 
|l5[ , and the spectrum of the Sinai billiard jl6| . 
The main conclusion is that if the corresponding 
classical system is chaotic, the microscopic cor- 
relations of the spectrum within an irreducible 
subspace are given by the RMT with the corre- 
sponding anti-unitary symmetry. 

4. CHIRAL RANDOM MATRIX MODEL 

In this section we will introduce a random ma- 
trix model for the QCD partition function. Our 
hope is to construct a model that reproduces both 
the microscopic spectral density and the correla- 
tions in the bulk of the spectrum. In the spirit 
of the invariant RMT we construct a model for 
the Dirac operator with the global symmetries of 
the QCD partition function as input, but other- 
wise gaussian random matrix elements. The chi- 
ral random matrix model (chRMM) that obeys 
these conditions is defined by p| , p7|Jl8| 

N f 

Z?= I DWjldetiV + m^e-^™^, (9) 
J f=i 

where 

and W is a rectangular n x m matrix with v = 
| n — rn | and N = n + m. The matrix elements 
of W are either real (/? = 1, chiral Gaussian Or- 
thogonal Ensemble (chGOE)), complex (j3 = 2, 
chiral Gaussian Unitary Ensemble (chGUE)), or 
quaternion real (/? = 4, chiral Gaussian Symplec- 
tic Ensemble (chGSE)). A precursor of this model 
was introduced within the framework of the in- 
stanton liquid partition function ]l9[] . This model 
reproduces the following symmetries of the QCD 
partition function: i) The £7,4(1) symmetry All 



nonzero eigenvalues of the random matrix Dirac 
operator occur in pairs ±A. ii) The topological 
structure of the QCD partition function. The ma- 
trix D has exactly \v\ = \n — m\ zero eigenvalues. 
This identifies v as the topological sector of the 
model. Hi) The flavor symmetry is the same as in 
QCD. For (3 = 2, it is SU(N f ) x SU(N f ). For /3 = 
1, it is SU(2N f ), and for @ = 4, it is SU(N f ). iv) 
The chiral symmetry is broken spontaneously for 
two or more flavors according to the pattern |20| 
SU(N f ) x SU(N f )/SU(N f ), SU(2N f )/Sp(N f ) 
and SU(Nf)/0(N f ) for /3 = 2, 1 and 4, respec- 
tively, the same as in QCD pl| . The chiral con- 
densate in this model follows from the Banks- 
Casher relation, S = limjv^oo Trp(0)/N. (N is in- 
terpreted as the (dimensionless) volume of space 
time.) v) The anti-unitary symmetries. For three 
and more colors with fundamental fermions the 
Dirac operator has no anti-unitary symmetries, 
and we choose = 2 in (§). For N c = 2 and 
fundamental fermions the Dirac operator satis- 
fies [CT2K, i'jD] = 0, where C is the charge con- 
jugation matrix and K is the complex conjuga- 
tion operator. Because, [Ct-iK) 2 = 1, the ma- 
trix elements of the Dirac operator can always be 
chosen real, and the corresponding random ma- 
trix ensemble is the chGOE. For two or more col- 
ors with fermions in the adjoint representation 
ijD has the symmetry [CK, ijD] = 0, but now 
(CK) 2 = —1, which allows us to rearrange the 
matrix elements of the Dirac operator into real 
quaternions. The corresponding random matrix 
ensemble is the chGSE. 

The ensemble of matrices (|l0]) weighted accord- 
ing to ([)]) is also known as the Laguerre ensemble. 
Note that its spectral correlations in the bulk of 
the spectrum are given by the invariant random 
matrix ensemble with the same value of /3 [p2| . 

5. UNIVERSAL PROPERTIES OF THE 
DIRAC OPERATOR 

5.1. Spectral correlations of the lattice 
QCD Dirac operator 

Recently, Kalkreuter calculated all eigen- 
values of the lattice Dirac operator both 
for Kogut-Susskind (KS) fermions and Wilson 
fermions for lattices as large as 12 4 . The ac- 
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Figure 1. Spectral correlations of Dirac eigenvalues for Wilson fermions (upper) and KS-fermions (lower). 



curacy of the eigenvalues was checked via sum 
rules for the sum of the squares of the eigenval- 
ues the lattice QCD Dirac operator. The anti- 
unitary symmetry of the KS and Wilson Dirac 
operator is different. For KS fermions we have 
@, [t 2 K, D ks ] = 0, with {t 2 K) 2 = -1, whereas 
for the Hermitean Wilson Dirac operator, 
[75C^T 2 ,7 57 C wilson ] = 0, with ( 75 CXr 2 ) 2 = 1. 
Therefore, we expect that the eigenvalue correla- 
tions are described by the GSE and the GOE, 
respectively [p5| . In Fig. |l| we show the re- 
sult for P(S), S 2 (n) and A 3 (n). The results 
for KS fermions are for 4 dynamical flavors with 
ma = 0.05 on a 12 4 lattice. The results for Wilson 
fermion were obtained for two dynamical flavors 
on a 8 3 x 12 lattice. The values of j3 and k are 
given in the label of the figure. For a discussion 
of other statistics we refer to |Q . 

An interesting question is the fate of spectral 
correlations for KS fermions in the continuum 
limit. To answer this question we have analyzed 
the 100-200 eigenvalues closest to zero. Even for 
the weakest coupling that was studied (/3 = 2.8) 
no deviation from the GSE was seen. 



5.2. Microscopic spectral density for a liq- 
uid of instantons. 

The investigation of the microscopic spectral 
density requires a very large number of indepen- 
dent gauge field configurations, which is difficult 
to generate by a lattice QCD simulations (see 
however |Ti| ). However, it is possible to simulate 
a large number of independent instanton liquid 
configurations. If the microscopic spectral den- 
sity of such 'smoothened' field configurations is 
given by chRMM, one certainly expects that this 
is the case for lattice QCD configurations. In Fig. 
|^ we show the spectral density of the Dirac oper- 
ator for N c = 2 and 3, and for 0, 1 and 2 massless 
flavors ]Ti| . In all cases do we find good agree- 
ment with the chRMM prediction. For N c = 3, 
Nf flavors and topological charge v it is given by 



PS(U) = - (J 2 a (u) - J a+ l(u)J a -l(u)) , (11) 

where a — Nf + v. The result for N c — 2, which 
is more complicated, is given in [ p8f . Qualita- 
tively, the microscopic spectral density has been 
observed for SU(2) staggered fermions and 
the lattice Schwinger model pj|. 
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Figure 2. The microscopic spectral density for a 
liquid of instantons. 

5.3. Valence quark mass dependence of the 
chiral condensate 

An alternative way to probe the Dirac spec- 
trum was introduced by the Columbia group 
p(i[ . They studied the valence quark mass de- 
pendence of the Dirac operator, i.e. £(m) = 
jr J d\p(\)2m/(\ 2 + to 2 ), for a fixed sea quark 
mass. In the regime (|l|), the valence quark mass 
dependence can be obtained analytically from the 
microscopic spectral density of M Q 



£(*) 



x(I a (x)K a (x) + I a+1 (x)K a ^(x)), (12) 



where x = mVE is the rescaled mass and a = 
Nf + v. In Fig. H we plot this ratio as a function 
of x for lattice data of two dynamical flavors with 
mass ma — 0.01 and N c = 3 on a 16 3 x 4 lattice. 
We observe that the lattice data for different val- 
ues of (3 fall on a single curve. Moreover, in the 
mesoscopic domain (|l|) this curve coincides with 
the random matrix prediction for Nf = v = 0. 
Apparently, the zero modes are completely mixed 
with the much larger number of nonzero modes. 
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Figure 3. The valence quark mass dependence of 
the chiral condensate. 

For eigenvalues much smaller than the sea quark 
mass, we expect to see the Nf = eigenvalue 
correlations. 

6. SCHEMATIC MODEL THE QCD 
PARTITION FUNCTION 

In this section we review results for a chiral 
random matrix model with nonzero temperature, 
T, and chemical potential, fi. The determinis- 
tic parts of this model are given by the matrix 
elements of jodo + njo in a basis exp[(i((2n + 
l)ir — &xgP)tT)](f) n {x) , where P is the Polyakov 
loop (introduced in Qj). If the matrix elements 
of the remaining terms are replaced by the ran- 
dom matrix ( |l0| ) we arrive at the model fl32|-]3l| 



V = 







iW + in T - 




(13) 



where £It — T <g)„ ((2n + 1)tt — argP)l. Inspired 
by |37f| , the simplest model is obtained by keeping 
only the lowest Matsubara frequency |32|33|. The 
random matrix model ( |l3| ) shows a phase transi- 
tion in T and fi. A related RMT as a schematic 
model of the Wilson Dirac operator at T = /x = 
was studied in ]3q] . Alternative random matrix 
models for chiral phase transition have been con- 
sidered as well. I only mention the study of the 
Dirac operator in the Gross- Witten model, where 
the transition is driven by the gauge field dynam- 
ics||. 

6.1. T^0 and fi = 

The spectral density of the model with only 
Matsubara frequencies ±nT can be obtained an- 
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alytically from the solution of a cubic equation 
It shows a second order phase transi- 
tion with critical temperature T c = 1 /ttT, (in the 
model with one Matsubara frequency) . Such 
a transition is typical for a 4-fermion model that 
is obtained after integrating over the random ma- 
trices. A particular interesting feature that could 
be checked by recent lattice studies of the QCD 
Dirac spectrum |27],[4(| , concerns the fluctuations 
of the smallest eigenvalue which is qualitatively 
different below and above T c p2| . 

If the phase of the Polyakov line is non-zero, the 
smallest Matsubara frequency is lower resulting 
in a higher critical temperature Q|. This result, 
that can also be obtained from an analysis of the 
NJL model @,^] , explains recent lattice data by 
the Columbia group [|30f showing that chiral sym- 
metry is restored later for the field configurations 
with a nonzero Z3 phase. 

Finally, the model (|l^) provides additional evi- 
dence for the universality of the microscopic spec- 
tral density (43). It can be shown that in 
spite of the dramatic change of the average spec- 
tral density, the microscopic spectral density does 
not change below T c . 

6.2. T = and /i ^ 

In this section we study the model (|l^) in 
the limit T -> at finite N. Note that the 
T — > limit taken after N — > oo is more 
subtle, e.g., the Fermi-Dirac distribution is ob- 
tained only after summing over all Matsubara 
frequencies EEf . For nonzero fi the Dirac oper- 
ator is no longer Hermitean, and its eigenvalues 
are scattered in the complex plane. An impor- 
tant question that can be asked in this context is 
the nature of the quenched limit. In agreement 
with earlier numerical results by Gocksch pc| , 
Stephanov |}5) showed analytically for the model 
( |l3| ) that the quenched limit is obtained as the 
limit limjVf-yO |det2?| JV/ . The absolute value can 
be interpreted as the limit of a model with both 
quarks and conjugate quarks which can produce 
a Goldstone boson with nonzero baryon number. 
This explains that fx c ~ \frn ||7j in quenched 
lattice QCD simulations p5[ . For more details, 
including results for the spectral density, we refer 
to Esll . A confirmation of some of these results 



can be found in (49j The failure of the quenched 
approximation at fi ^ was also observed for a 
one-dimensional U(l) model J50|] . 

In the unquenched case, the random matrix 
model shows a phase transition at /i c = 0.53 
|35|,[36|]. By putting the determinant inside the 
operator, the resolvent of the Dirac operator, 
G(z) — Tr(z — D)~ 1 , can be obtained numeri- 
cally. For z inside the domain of the eigenvalues, 
it diverges in the thermodynamic limit, and is dif- 
ferent from the quenched result obtained analyti- 
cally in p5fl . For z outside this domain quenched 
and unquenched results coincide in the thermo- 



dynamic limit |36|. The same phenomenon is ob- 
served H for the U(l) model of @. 

7. DISCUSSION AND CONCLUSIONS 

We have shown that the spectrum of the QCD 
Dirac operator shows universal features that can 
be obtained from a random matrix model with 
the global symmetries of QCD. In this way, we 
have obtained analytical results for the finite vol- 
ume corrections to the valence quark mass de- 
pendence of the chiral condensate and the spec- 
tral correlations in the bulk of the spectrum. We 
have also shown that an extension of this random 
matrix model provides a schematic model for the 
chiral phase transition. Interesting results have 
been obtained for QCD at nonzero temperature 
and at nonzero chemical potential. 

This work was partially supported by the US 
DOE grant DE-FG-88ER40388. 
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